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Abstract 

The buckling characteristics of Reng 41 tubu- 
lar panels for a hypersonic aircraft wing were 
investigated. The panels were repeatedly tested 
for buckling characteristics using a hypersonic 
wing test structure and a universal tension/ 
compression testing machine. The nondestructive 
buckling tests were carried out under different 
combined load conditions and in different tem- 
perature environments. The force/stiffness tech- 
nique was used to determine the buckling loads of 
the panels. In spite of some data scattering 
resulting from large extrapolations of the data- 
fitting curve (because of the termination of 
applied loads at relatively low percentages of the 
buckling loads), the overall test data correlate 
fairly well with theoretically predicted buckling 
interaction curves. Also, the structural effi- 
ciency of the tubular panels was found to be 
slightly higher than that of beaded panels. 


f cr stress intensity at failure 

f C y compression yield stress 

fpi proportional-limit stress 

(threshold of inelastic 
stress region) 

G shear modulus 

1 index 



k F(703) cr /F ( 703)*, an extrapolation 

factor 

m exponent in the expression of D 


D 

Dcr 


Nomenclature 

Nx 

panel axial compression load, lb/in 

coefficients, where j = 1 to 6 

(N x )cr 

kN x , the value of N x at buckling. 

generalized strain variable 


lb/in 

buckling-critical generalized strain 
variable 

★ 

N X 

the value of N x at maximum applied 
load 


E 

E sec 

Etan 

F 


modulus of elasticity 
secant modulus 
tangent modulus 
applied load 


N X y panel shear load, lb/in 

(^xy)cr KN X y, the value of N X y at buckling, 

lb/in 

* 

N X y the value of N X y at maximum applied 

1 oad 


F* maximum applied load 

F cr critical load 

F/S force/stiffness 

F ( 7 03 ) applied force at HWTS location 703 


n 

P 

R 


shape factor in Ramberg-Osgood 
stress-strain law 

lateral pressure 

radius of circular arc of the panel 
tube cross section 


F(703) cr kF ( 703)* , the value of F( 703) at 
predicted buckling point 


F(703)* maximum value of F(703) in the 

nondestructive buckling test 


fcb> fcc> 
fcs 


bending, compression, and shear local 
buckling-critical stresses of 
circular elements 


RSG-j 

RSG? 

RSG-, 


rosette strain gage output 


rosette strain gage output at 
maximum applied load F* 

rosette strain gage output with 
structure at uniform elevated 
temperature without applied 
mechanical load 
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RSGijk rosette strain gage ljk, or its 

output 

S c 2aR , the arc length of circular arc 

element of the panel tube cross 
section 



SGi axial strain gage output 


SGt axial strain gage output with 

structure at uniform elevated 
temperature without applied 
mechanical load 

T temperature 

t thickness of the tubular wall 


t effective thickness of the panel 

w unit panel weight 

S c 

a — -, one-half the angle of circular 

c K 

arc of the panel tube 
y shear strain 

Ter critical buckling strain in shear 

eb bending strain 

( e b) cr critical buckling strain in bending 

ec axial compression strain 

( e c) cr critical buckling strain in 

compression 

n plasticity correction factor, where 

1 1 = 1 to 3 

o tan plasticity correction factor 

a calculated from tangent 

modulus 

n plasticity fcorrection factor 

calculated from secant 
modulus 

Oj itan> the plasticity correction 

factor calculated from 
tangent modulus 


structural temperature to 350°F (a warm structure). 
Another concept proposed for future hypersonic 
aircraft* is an aerodynamically acceptable wavy 
heat shield made of heat-resistant metal (such 
as Rene 41) designed to limit the structural tem- 
perature to about 1350°F (a hot structure). 

Buckling is a major concern for hot struc- 
tures because of the combined effects of thermal 
stresses and air loads. For future hypersonic 
flight vehicles, advanced structural concepts that 
show promise of low structural unit mass and high 
buckling strength have been identified. 2- ** 

Because curved shell sections exhibit high local 
buckling strength, most of the structural panel 
concepts investigated apply curved surfaces to 
achieve high buckling strength. Two of the hot 
structural panel concepts are beaded panels and 
circular tubular panels. Results from extensive 
buckling studies for Ren£ 41 beaded panels and / 
aluminum circular tubular panels are reported in 
Ref. 9 and Refs. 2 to 6, respectively. However, 
all the results for tubular panels were obtained 
from machine-controlled, room-temperature loading 
tests, not from tests of panels in a heated, 
built-up wing structure. Furthermore, only 
limited buckling data have been reported on 
Ren§ 41 noncircular tubular panels, which, for 
lightly loaded conditiqns, are more efficient than 
circular tubular panels, based on the computer 
codes analysis methods. 6 

Thus, to characterize the buckling behavior 
of the tubular panels in a heated wing structure, 
five Ren£ 41 noncircular tubular panels were 
attached to the wing root region of the hypersonic 
wing test structure (HWTS, described in the Test 
Equipment section) for extensive nondestructive 
buckling tests under different combined load con- 
ditions (axial compression, bending under lateral 
pressure, and shear). The temperature environ- 
ments used in the tests were 70°F, 550°F, and 
1000°F. The buckling loads were estimated by 
using the force/stiffness (F/S) method of plot- 
ting the test data, 10 and the results were com- 
pared with theoretically predicted buckling 
interaction curves. 


Test Panels 


r)^ i S ec» ttle plasticity correction 

factor calculated from 
secant modulus 

the plasticity correction 

factor for axial compression 
and bending 

v Poisson's ratio 

Introduction 

When an aircraft reaches the hypersonic speed 
range, aerodynamic heating becomes significant. 
Hypersonic flight vehicle concepts that have been 
advanced include either using a thermal protection 
system (TPS) to prevent the vehicle from over- 
heating or allowing an exposed vehicle surface 
to reach a radiation equilibrium temperature with 
respect to aerodynamic heating. For example, the 
space shuttle uses a TPS designed to limit the 


Each test panel was fabricated from two sheets 
of Ren§ 41 alloy seam-welded together to form five 
flat regions (double sheets) and four noncircular 
tubular regions (flattened tubes formed by the 
union of two circular arcs). The dimensions of 
the test panel and the geometry of the noncircular 
tube cross section are shown in Fig. 1. Figure 2 
is a photograph of one of the six tubular panels 
that were fabricated. The material was procured 
in a solution-annealed (1975°F) condition. Prior 
to welding, the circular arcs in each sheet were 
incrementally brake formed, and the end closures 
were die formed. Doublers were spot welded to the 
upper and lower ends of the panel to prevent local 
end failure and to reduce excessive deformations 
due to shear loads. After final weld assembly, 
the panels were aged for 1 fir at 1650°F followed 
by 10 hr at 1400°F. Each panel has eight attach- 
ment points for Z-shaped clips to support the 
heat shields, which will be described later. 

A detailed description of the fabrication pro- 
cess is given in Refs. 3 and 6. 
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Analysis 

Local Buckling 

For a tubular panel with tubes of complete 
circular cross section, the equations for local 
buckling (bead crippling) of circular arc elements 
of the panel in compression, bending, and shear 
may be written as follows 2 

For compression (Ref. 2, eq. 14-3), 

f cc = O.738n 3 E0 1-19 (1) 

For bending (Ref. 2, eq. 12-33), 

f C b = O.77n 3 E0 K15 (2) 

For shear (Ref. 2, eq. 12-34), 

f cs = n 2 GK s^|~^ ( 3 ) 


Eqs. (2), (3), and (5), are defined as (Ref. 2, 
p. 31) 


n l 

= n tan = ^tan/E 

(6) 

n 2 

= n sec = ^sec/E 

(7) 

n 3 

= n l n 2 = n tan n sec 

(8) 

E tan 

fcr 


f cr/E + n[0.002 (f cr /f C y) n 

- 0.00001] 

E sec 

_ fcr 

(9) 

f cr /E + 0.002 (f C r/fcy) n - 

0.00001 



(10) 

where 

fcr > fpt = f C y ( 0 . 005 ) 1/n 

(11) 


and 


where 


^tan - E sec - E 


if f cr < fpjj 


Ks 




(4) 


Buckling equations (1) and (2) are valid for the 
range 20 < R/t < 120, and Eq. (3) is valid for 


5^ 

— Vl - v 2 > 50 
Rt 


To apply Eq. (1) to the noncircular tubular panel, 
a knockdown factor of 0.86 is needed (see Ref. 6, 
p. 46). Equations (2) and (3) are applied 
directly to the flattened tubular panel without 
modification. The buckling equation for compres- 
sion of the noncircular tubular panel is then 


f CC = (O.86)(O.738)n 3 E0 1,19 

= 0.635n 3 E^y* 19 (5) 

Equations (2), (3), and (5) were used to calcu- 
late the theoretical buckling strains in bending, 
shear, and compression, respectively, for use in 
the F/S plots of the nondestructive buckling data 
described in the following section. The plasti- 
city correction factors n-, and n 3 , which appear in 


The shape factor n in tfie Ramberg-Osgood stress- 
strain laws for Reng 41 (Eqs. (9) to (11)) was 
taken as n = 25 at 70°F, 22.2 at 550°F, and 18.5 
at 1000°F. Figures 3, 4,. and 5 show the plots of 
equivalent elastic stress f C r/n 1 (i = 1,2,3) as a 

function of^actual stress f cr for the three tem- 
peratures 70°F, 550°F, and 1000°F, respectively. 
Figure 6 shows values of the modulus of elasticity 
E and shear modulus G for Ren£ 41 as functions of 

temperature. 12 With the aid of Figs. 3 to 6, the 
theoretical critical buckling strains in compres- 
sion (e c ) cr , m bending (cb)cr» and m shear y cr 
can be calculated from Eqs. (5), (2), and (3), 
respectively, as 

( £ c)cr 


( £ b)cr 


Ter 


fee 

n 3 E 

(12) 

feb 

n 3 E 

(13) 

^CS 

^2^ 

(14) 


The values of (e c ) cr , (et>) cr , y C r (with corre- 
sponding values of n 2 and n 3 ), E, and G for the 
different temperatures are tabulated below. 





70°F 




550 c 

'F 



1000°F 

( c c)c r 

2.330 

X 

10-3 

(n 3 = 1) 

2.330 

X 

10-3 

(n 3 = 1) 

2.330 

X 

10-3 ( n3 = i) 

( £ b )cr 

3.412 

X 

10-3 

(n 3 = 0.96) 

3.412 

X 

10-3 

(n 3 = 0.97) 

3.412 

X 

10-3 (n 3 = o.99) 

Ter 

2.208 

X 

10-3 

(n 2 = 1) 

2.208 

X 

10-3 

(n 2 = 1) 

2.208 

X 

10-3 ( n2 = i) 

E (1 b/i n2) 

31.6 

X 

106 


29.1 

X 

106 


26.5 

X 

106 

G (lb/in 2 ) 

12.1 

X 

10 6 


11.2 

X 

106 


10.2 

X 

106 
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Notice that only the bending buckling strains are 
in the plastic range. 

F orce/Stiff ness Method 

The purpose of conducting nondestructive 
buckling tests instead of destructive buckling 
tests was to avoid the cost associated with 
destructive tests of a large number of panels. 

In destructive buckling tests, only one buckling 
data point for a certain load condition can be 
generated from each test panel. However, by using 
the F/S method to predict the buckling strength, 
a wide range of buckling data points can be gen- 
erated from the same set of test panels. The F/S 
method described in this paper was advanced by 
Jones and GreeneH for the prediction of buckling 
strength of complex structural components having 
buckling behavior that is quite complex (for 
example, local buckling) or nonlinear. 

Although local instability is usually the 
critical failure mode for the tubular panels, 
local and panel-general instability may occur 
at similar loading levels when lateral pressure 
is low and shear is high (see Results section). 

If the panel -general instability were to occur 
first, local instability would probably also 
occur. Consequently, the F/S method used to 
identify buckling strength was one that predicts 
local buckling failure. The method uses a plot 
of force F against stiffness F/D, where F can be 
any particular load of the applied load set and 
D is a generalized strain variable that accounts 
for axial compression, bending, and shear com- 
ponents. The generalized strain variable D is 
given by 


the F/S method depends on the following factors' 
location of strain gages (they must measure 
strain that is sensitive to the impending buck- 
ling mode shape), distance of extrapolation 
(that is, how far the test data stopping point 
falls from the limit strain line), how well the 
curve fits the data, and the accuracy with which 
the limit strain line itself is determined (for 
example, critical strains determined empirically 
for a specific configuration may be more accurate 
than critical strains determined analytically from 
general equations). 

Extrapolation of the F/S data-fitting curve to 
its intersection with the limit strain line is 
accomplished through least squares fitting of the 
test data by using the following third-order 
equation 

F = 1 + C t F + C 2 f2 ( 17) 

0 " C 3 + C 4 F + C 5 f2 + C 6 F 3 

Thp critical value of F (the intersection 
point F cr ) is determined by setting D = 1 in 
Eq. (17). In this F/S analysis, D is expressed 
as follows 

For room temperature, . 

_ SGi + SG 2 + | SGj - SG 2 | + 2 | RSG 2 - RSG 3 | 2 
3 ( £ c)cr 2(£b)cr V3ycr 

(18) 

For elevated temperatures, 


D 


e c + £ b 
( £ c)cr ( e b)cr 


+ 



(15) 


(SGj - "5G"i) + (SG 2 - 3U 2 ) 
2( £ c)cr 


and the critical local buckling condition occurs 
when 


D - Dcr = 1 


(16a) 


and 


+ K SG 1 - SG l) ~ ( SG 2 - SG 2) 

2( e b)cr 


2 | (RSG 2 - RSG 2 ) - (RSG 3 - RSG 3 ) I 2 
V3y cr 


The strains e c , e(d, and y are measured with strain 
gages, and the critical strains (e c ) cr , ( £ b)cr’ 
and y cr (tabulated previously) are calculated from 
Eqs. (12), (13), and (14), respectively. The 
exponent m in Eq. (15) was empirically determined 
to be 2 for most types of panels, including the 
complete circular tubular panels. 5 For the F/S 
analysis described here, the value m = 2 was used. 


(19) 

The values SGj and SG 2 are the outputs of the 
axial strain gages placed, respectively, at the 
lower and upper outermost fibers of the tube at 
the panel center region and RSG 2 and RSG 3 are the 
outputs of the rosette strain gage legs oblique 
to the tube axis. Finally, the overbar indicates 
the initial nonzero strain gage readings at ele- 
vated temperatures when no mechanical loads were 
applied (panels were soaked at uniform temperature 
to determine these initial strains). 


Equations (15) and (16) represent a critical 
strain interaction surface that is the basis for 
the limit strain lines used in the F/S plots. 
Figure 7 graphically illustrates the F/S method, 
which requires extrapolation of the curve fitting 
the test data points. The buckling failure load 
is determined from the intersection of the extrap- 
olated curve and the limit strain line. The 
accuracy of buckling failure load prediction using 


Panel B uc kling Loads 

After the critical load F cr is determined 
using the F/S method, the associated critical 
panel axial compression load (N x ) cr and critical 
panel shear load (N X y) cr must be determined. If 
Nj and N X y are the panel axial and shear loads, 
respectively, associated with the maximum applied 
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load F* (see Fig. 7) and if (RSG* - RSG 1 ) , 

(RSG* - RSG 2 ), and (RSG^ - RSGj) are the readings 
of the three legs of the rosette strain gage 

when F = F*, ((RSG* - RSGj) being in the axial 
★ ★ 

direction), then N x and N X y can be calculated as 
N x = Et(RSG* - RSGj) (20) 

N X y = ~ Gt(2t ) | (RSG 2 - RSG 2 ) - (RSG* - RSG 3 ) | 

( 21 ) 

where t is the effective thickness of the panel 
and RCG* (l = 1,2,3) are the rosette strain gage 
readings at F = F*. 

If the extrapolation factor k (see Fig. 7) is 
defined as 



then (N x ) cr and ( N X y ) c r can be estimated as 

(N x ) C r = kN* (23) 

( N xy )cr = ^N X y (24) 

The values of (N x ) cr and (N X y) cr thus obtained are 
used in constructing the buckling interaction 
plots. 

Test Equipme nt 

Combined L oads 

Hypersonic Ming T est Structure. The hyper- 
sonic wing test structure (HWTS, 85 ft2 in area) 
shown in Fig. 8 is a portion of the wing of a pro- 
posed hypersonic research airplane (HRA). shown in 
Fig. 9. The HRA is to be built for cruising at 
Mach 8 with 2.5 g pullup capability (Ref. 13 
discusses the HRA in great detail). The HWTS was 
constructed based on the knowledge gained from the 
hot structural concepts investigations in the 
hypersonic cruise vehicle research program. 1 * 14 
The HWTS has been tested extensively to evaluate 
the structural concept itself, the final struc- 
tural design, the flight loads instrumentation, 
high-temperature calibration methods, and tem- 
perature simulation techniques. The beaded skin 
panels and corrugated spars and ribs are made of 
Rene 41 nickel alloy. The heat shields are also 
made of Ren£ 41 alloy except for those along the 
leading edge, which are made of TD Ni-20Cr. The 
heat shields are single-sheet panels, which are 
slightly corrugated in the chordwise direction and 
capable of operating with surface temperature in 


excess of 1800°F. The heat shields are separated 
from the beaded skin panels by Z-shaped support 
clips to minimize heat conduction from the heat 
shields to the substructure. The HWTS is mounted 
inverted to the support structure through a tran- 
sition section (Fig. 8), so the compressive sur- 
face of the actual vehicle would be the lower 
surface of the HWTS. The transition section is 
included to provide a buffer between the support 
structure and the test portion of the wing. The 
lower wing root zone is the most critically 
compression-loaded area, and the five beaded 
panels there can be replaced with different types 
of panels for testing structural efficiencies. 

The five tubular panels were installed in this 
lower wing root area for the nondestructive 
buckling tests. Figure 10 shows the HWTS with the 
heat shields removed to reveal the substructure 
and the Z-shaped clips that support the heat 
shields. 

Mechanica l Loadin g Systems . Figure 11 shows 
the locations of the five test tubular panels and 
the applied mechanical load distribution on the 
HWTS. Twenty channels of closed-loop electro- 
hydraulic equipment were used to apply mechanical 
loads to the test structure at different load 
points, as shown in Fig. 11. Ten hydraulic jacks 
were used to apply vertical loads (simulating lift 
load) to the HWTS to irfduce compression loads to 
the test panels, eight of these jacks applied 
loads through two-point whiffletrees. Horizontal 
loads (simulating drag and thermal load from fore 
and aft wing structure) were applied by the 
remaining 10 hydraulic jacks at single load points 
on the fore and aft edges of the HWTS to induce 
shear loads to the test panels. The pressure 
loads (which induced bending loads) normal to the 
upper surfaces of the test panels were generated 
through the stainless steel pressure pans posi- 
tioned over the test panels. 

Heating System . For the HWTS elevated- 
temperature tests (simulating aerodynamic 
heating), a specially constructed heating system 
was used (Fig. 12). The system was designed to 
heat the entire upper and lower surfaces of the 
HWTS to temperatures corresponding to a Mach 8 
flight profile. Infrared quartz lamps mounted on 
polished aluminum reflectors (Fig. 12) were used 
to provide the radiant heat flux for the simula- 
tion of aerodynamic heating. The reflectors were 
water cooled to enable them to withstand the high 
temperatures achieved during the heating tests. 
Both the lower and the upper heating systems were 
slightly contoured to match the surface shape of 
the HWTS and were mounted on rollers and tracks 
(see Fig. 12) so they could be easily removed for 
access to the HWTS and then precisely reposi- 
tioned. The heating systems were positioned with 
the reflector surfaces approximately 6 in. from 
the heat shields of the HWTS. Because the test 
structure has mechanical load attachment points on 
the lower surface, gaps in the lower heater were 
provided along the spar caps to enable con- 
nections with the hydraulic vertical loading 
system (Fig. 13), these gaps between load points 
were filled by a double row of quartz lamps 
mounted on separate, long, narrow, water-cooled 
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aluminum reflectors (strip heaters) installed 
parallel to the spar caps. Heating rates were 
controlled through thermocouples attached to the 
test structure exterior surfaces. During the 
heating tests, insulation curtains (see Fig. 12) 
could be draped around the HWTS and the heating 
systems to reduce radiation and convection heat 
losses. The plumbing for the reflector cooling 
water included a pressure gauge for each feed line 
to ascertain cooling provisions. 

Compr e ssion Only 

For gathering additional room-temperature 
(70°F) buckling data in pure compression, a uni- 
versal tension/compression testing machine was 
used in the individual panel compression buckling 
tests. Figure 13 shows the test machine with the 
test panel mounted. The surfaces of the upper and 
lower platens (which come into direct contact with 
the panel edges) were machined flat to ensure pure 
compression loading and eliminate possible bending 
effect due to misaligned platen surfaces. The 
lower platen was rested on a spherical seat and 
provided proper alignment with the test panel. 

End supports mounted on the panel provided 
surfaces for load transfer and served as rein- 
forcement to eliminate warping of the panel ends. 
The surfaces of both end supports were milled 
parallel to each other and perpendicular to the 
panel tube axes to provide pure compression load 
transfer. 

The panel vertical edges were bolted to the 
Z-section stiffeners to approximate the edge 
stiffness conditions of the panel if it were 
mounted on the wing. The interfaces between the 
panel and the Z stiffeners were lubricated. The 
bolt holes on the Z stiffeners were oval in shape 
so the bolts could move when the panel deformed. 

A total of 11 displacement transducers (DT) were 
used to measure the out-of-pl'ane deformations of 
the test panel at 11 locations (Fig. 13). 

Strain Gage Instrumentation 


acquisition system for strain gage measurements 
was ±4.88 microstrain, which represents 0.3 per- 
cent of the strain gage calibrated output. This 
provided accurate recorded data. Figure 15 shows 
the fully instrumented test panels attached to 
the HWTS lower wing root test area with panel 3 
removed to show the pressure pan interior. 

Tests 

Combined Load Condi t ions 

To generate a wide range of buckling data for 
buckling interaction plots using the F/S method, 
a series of nondestructive buckling tests was 
conducted under various combined load conditions 
at three temperature levels (70°F, 550°F, and 
1000°F). Table 1 shows the maximum loads applied 
at the load points for different load conditions. 
The pressure load was always maintained at the 
constant level of 0.75 1 b/i n 2 or at zero. For 
load conditions 4.1, 7.1, and 8.1, which do not 
have horizontal loads, the test panels were under 
very little shear load. Before the series of 
tests at room temperature (70°F), a pressure 
check was conducted to ensure that a constant 
pressure level could be maintained during the 
tests. For the elevated-temperature tests, a 
heating check was conducted before the actual 
buckling test to ensure- that constant temperature 
levels could be maintained over one test period. 
Table 2 shows test identification numbers and 
corresponding load conditions. 

Compress ion-Only Tests 

Both panels 1 and 3 were used for single-panel 
room-temperature (70°F) compression tests up to 
buckling failure. During the tests, the outputs 
from the load cell, strain gage channels, and DT 
channels were fed into the data acquisition system 
so that F/S plots could be generated. The buck- 
ling loads obtained from the F/S plots could then 
be compared with the actual buckling loads. 

Data Reductions 


The strain gage locations on the surfaces of 
the five test panels are shown in Fig. 14, looking 
downward onto the top of the test panels. The 
strain gages with ( ) symbols were located on 
the upper surfaces of the panels, the rosette 
strain gages with [ ] symbols were used for the 
elevated-temperature tests. The strain gages on 
the upper and lower surfaces of the panel tubes 
are axial strain gages of two types foil type 
(circular symbol) and capacitance type (square 
symbol). Two of these strain gages, one of each 
type, were mounted end to end at each location 
shown in Fig. 14. Of the axial strain gages, 
only the capacitance strain gages were capable of 
operating at temperatures above 550°F. The strain 
gages on the surface of the panel flat area were 
of the delta-rosette foil type (indicated by 
triangles in Fig. 14). The delta-rosette strain 
gages were used to make measurements at three 
angular orientations, spaced 120° apart, starting 
in the direction parallel to the wing spars and 
rotating clockwise (when looking down on the test 
panels and inboard). The accuracy of the data 


In applying the F/S method described pre- 
viously, a typical vertical load (lift force), 
F(703) (located at load point 703, associated 
with jack number 3, Fig. 11) was chosen for F in 
Eq. (16) (or Fig. 7) in the F/S calculations for 
all the test panels. The force F(703) was arbi- 
trarily selected as a representative measurement 
of all applied loads since all loads were directly 
proportional to each other and were applied simul- 
taneously. For F = F ( 703 ) , Eq. (22) becomes 


Fcr = F ( 703 ) c r 
F* F (703)* 


(25) 


where F(703) cr is the value of F ( 703 ) at the pre- 
dicted buckling point and F ( 703 )* is the maximum 
value of F ( 703 ) in the nondestructive buckling 
test. The critical values of the panel loads 
(N x ) cr and (N X y) cr may be calculated by using 

Eqs . (20), (21), (23), and (24). For example, 
for panel 1 at 70°F with the rosette strain gage 
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combination RSG933, RSG934, and RSG935 (Fig. 14), 
Eqs. (20) and (21) may be written as 

N x = Et (RSG933* - 0) (26) 

Njy = ~ G(2t) | (RSG934* - 0) - (RSG935* - 0) | 


(27) 

and (N x ) cr and (N X y) cr can be calculated from (see 
Eqs. (23) and (24)) 

(Nx)cr = kN* = N* (28) 

(Nxy)cr = kN* y = N * y (29) 


Results 

Figures 16 and 17, respectively, show the F/S 
plots for two typical test cases 4.4.6 (p * 0) 
and 4.3.4 (p = 0). The solid curves shown in the 
figures were drawn by least squares fittings of 
the test data points. The least-squares-fit 
curves based on Eq. (17) tend to bend upward 
gradually far outside the range of the data 
points. For some test cases, or for certain 
strain gage combinations in the same test case, 
the least-squares-fit curves start to bend up- 
ward immediately after the last data points and 
intersect the limit strain lines at points that 
give rather high values of the buckling loads. 

For such cases, the extrapolations of the test 
data curves were accomplished by observation (or 
visual fittings). When there is lateral pressure 
(p * 0), the F/S curves are usually convex upward 
(Fig, 16). However, when there is no lateral 
pressure (p = 0), the F/S curves are usually 
strongly convex downward, except for the low- 
load region, giving quite accurate buckling load 
predictions (Fig. 17). The F/S plots for the rest 
of the test cases are similar and therefore are 
not shown. 

Figure 18 shows the F/S plots for panel 3 
under the room-temperature (70°F) single-panel 
compression tests. Note that the actual buck- 
ling points fell in the vicinity of the limit 
strain lines. This indicates that the F/S method 
could give relatively accurate predictions of 
buckling loads. 

In most of the F/S plots, the extrapolation 
factor k was in the approximate range of 2 to 3. 

Figures 19 to 23 show the buckling interaction 
plots for different temperature environments, with 
or without pressure. In the figures, theoretical 
buckling interaction curves (solid curves for 
panel instability, dashed curves for local buck- 
ling) are shown for comparison. These theoret- 
ical buckling interaction curves were generated 
through temperature and material modulus correc- 
tions of the results given in Ref. 6 (table 1 
(modified analysis), p. 48) which were calculated 
by using equations given in Ref. 2. The maximum 


applied load was limited to the values shown in 
Figs. 19 to 23 to prevent spar cap failure in the 
HWTS itself. Thus, large F/S data curve extrapo- 
lations were necessary. For high-compression 
loadings (Fig. 20), the maximum applied load 
levels were not high enough to give more accurate 
buckling data through the F/S data extrapolations. 
Nevertheless, the correlations between the test 
data and the predictions are fairly good in spite 
of data scatter due to the large extrapolations 
of the data-fitting curve for some test cases 
(for example, near the N x axis for T = 70°F, 
p = 0.75 lb/in 2 , Fig. 20). Most of the data 
points, including the actual buckling points 
obtained from the single-panel compression tests 
(Fig. 19), fall outside the predicted interaction 
curves (surfaces), indicating that the theory is 
conservative. In all other cases, the theoretical 
curves fall within the scatter of the experimental 
data. Therefore, the experimental buckling data 
verify the theory for the applied test conditions. 
Notice that the existence of lateral pressure 
decreases the compressional buckling load (N x ) cr 
considerably but decreases the shear buckling 
load (N X y) cr only slightly. Also, the buckling 
interaction curve shrinks as the temperature is 
increased because of decreases in E and G. 

Finally, the buckling strength per unit panel 
weight, (N x ) cr /w (w = 0.0161 lb/in 2 ), for 
panels 1 and 3 were (N x ) cr /w = 1.5068 x 10$ in 
and (N x ) cr /w = 1.3280 x 10$ in, respectively. 

These values are slightly higher than 
(N x ) cr /w = 1-1507 x 105 m (w = 0.0146 lb/in 2 ) 
for the beaded panels reported in Ref. 9. 

Conclusions 


Five Ren§ 41 tubular panels were simulta- 
neously and repeatedly tested for buckling charac- 
teristics using the hypersonic wing test 
structure. Nondestructive buckling tests were 
carried out under different combined load con- 
ditions and in different temperature environments. 
The force/sti ffness technique was used to deter- 
mine buckling loads for a wide range of load con- 
ditions without failing the test panels. In spite 
of some data scattering due to large extrapola- 
tions of the data-fitting curve, the overall test 
data correlate fairly well with theoretically pre- 
dicted buckling interaction curves. The existence 
of lateral pressure decreased the compressional 
buckling load (N x ) cr considerably. However, the 
effect of lateral pressure on the reduction of the 
shear buckling load (N X y) cr was quite small. 

Also, increasing the temperature was found to 
decrease both (N x ) cr and (N X y) cr because of reduc- 
tions in the values of E and G at elevated tem- 
peratures. The fact that almost all the test data 
for nearly pure compression (N X y small) at room 
temperature (70°F) fell outside the predicted 
buckling interaction curves indicates that the 
theory used to design and predict the buckling of 
panels is conservative for that condition. Other- 
wise, the test data verify the theory for all 
other test conditions. Finally, the tubular 
panels were found to be slightly more structurally 
efficient than beaded panels. 
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Table 1 Maximum loads applied at load points for different load combinations 


(a) T = 70°F 


Jack 

number 

Load cell 
size, lb 

Jack 

position 


Load 

condition 

identification 

numbers 


3.1 b 

3.6 C 

3.8 C 

4. ic 

4.2C 

4.3C 

4.4C 



Maximum load.d lb 



1 

5 x io3 

Vertical 

1680 

2521 

4,201 

3500 

2800 

1200 

-4000 

2 

10 

Vertical 

869 

1304 

2,173 

6000 

2800 

2000 

2400 

3 

10 

Vertical 

624 

935 

1,559 

6000 

2800 

2000 

2400 

4 

5 

Vertical 

863 

1294 

2,157 

3500 

2800 

-2800 

-4000 

5 

5 

Hon zontal 

0 

1650 

2,750 

0 

-5000 

-5000 

-5000 

6 

10 

Hori zontal 

0 

4260 

7,100 

0 

-6500 

-6500 

-6500 

7 

10 

Horizontal 

0 

-4030 

-6,716 

0 

6500 

6500 

6500 

8 

10 

Hori zontal 

0 

-1232 

-2,054 

0 

6500 

6500 

6500 

9 

10 

Hon zontal 

0 

-3840 

-6,350 

0 

6000 

6000 

6000 

10 

10 

Horizontal 

0 

737 

1,229 

0 

6000 

6000 

6000 

11 

10 

Vertical 

2265 

3398 

5,663 

3000 

800 

-800 

-4000 

12 

5 

Vertical 

-95 

-143 

-238 

1800 

1400 

1200 

1200 

13 

20 

Hon zontal 

0 

-6692 

-11,154 

0 

6500 

6500 

6500 

14 

10 

Horizontal 

0 

-3945 

-6,575 

0 

6500 

6500 

6500 

15 

10 

Horizontal 

0 

-1740 

-2,900 

0 

-6500 

-6500 

-6500 

16 

10 

Horizontal 

0 

2400 

4,000 

0 

-6500 

-6500 

-6500 

17 

5 

Vertical 

1267 

1901 

3,168 

3500 

2800 

2800 

4000 

18 

5 

Vertical 

-1846 

-2770 

-4,616 

3500 

4800 

4800 

4000 

19 

5 

Vertical 

1760 

2640 

4,400 

2800 

2800 

2800 

4000 

20 

5 

Verti cal 

542 

814 

1,356 

1800 

1400 

1200 

1200 


a Load number 0.75 was a 0.75-lb/in2 pressure test. 

^Pressure load maintained at 0. 

c Tests conducted at pressure loads of 0 and 0.75 lb/in2. 

^Positive values indicate tension, negative values indicate compression. 
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Table 1 Continued 


(b) T = 550°F 


Jack 

number 

Load cell 
size, lb 

Jack 

position 

Load condition identification number a > b 

6.6 

6.8 

7.1 

7.2 

7.3 

7.4 


Maximum load, c lb 



1 

5 x 103 

Vertical 

2521 

4,201 

3500 

2800 

1200 

-4000 

2 

10 

Vertical 

1304 

2,173 

6000 

2800 

2000 

2400 

3 

10 

Vertical 

935 

1,559 

6000 

2800 

2000 

2400 

4 

5 

Vertical 

1294 

2,157 

3500 

2800 

-2800 

-4000 

5 

5 

Horizontal 

1650 

2,750 

0 

-5000 

-5000 

-5000 

6 

10 

Hori zontal 

4260 

7,100 

0 

-6500 

-6500 

-6500 

7 

10 

Horizontal 

-4030 

-6,716 

0 

6500 

6500 

6500 

8 

10 

Horizontal 

-1232 

-2,054 

0 

6500 

6500 

6500 

9 

10 

Horizontal 

-3840 

-6,350 

0 

6000 

6000 

6000 

10 

10 

Horizontal 

737 

1,229 

0 

6000 

6000 

6000 

11 

10 

Vertical 

3398 

5,663 

3000 

800 • 

-800 

-4000 

12 

5 

Vertical 

-143 

-238 

1800 

1400 

1200 

1200 

13 

20 

Horizontal 

-6692 

-11,154 

0 

6500 

6500 

6500 

14 

10 

Hon zontal 

-3945 

6,575 

0 

6500 

6500 

6500 

15 

10 

Horizontal 

-1740 

-2,900 

0 

-6500 

-6500 

-6500 

16 

10 

Horizontal 

2400 

4,000 

0 

-6500 

-6500 

-6500 

17 

5 

Vertical 

1901 

3,168 

3500 

2800 

2800 

4000 

18 

5 

Vertical 

-2770 

-4,616 

3500 

4800 

4800 

4000 

19 

5 

Vertical 

2640 

4,400 

2800 

2800 

2800 

4000 

20 

5 

Vertical 

814 

1,356 

1800 

1400 

1200 

1200 


a Load number 550 was a 550°F heating check. 


Pressure 
c Positi ve 


load maintained at 0.75-lb/in2. 

values indicate tension, negative values indicate compression. 
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Table 1 Concluded 


(c) T = 1000°F 


Jack 

number 

Load cell 
size, lb 

Jack 

position 

Load condition identification number 3 

9.3b 

9.5b 

8.ic 

8.2C 

8.3C 

8 . 4 c 


Maximum load, d lb 



1 

5 x io3 

Vertical 

2521 

4,201 

3500 

2800 

1200 

-4000 

2 

10 

Vertical 

1304 

2,173 

6000 

2800 

2000 

2400 

3 

10 

Vertical 

935 

1,559 

6000 

2800 

2000 

2400 

4 

5 

Vertical 

1294 

2,157 

3500 

2800 

-2800 

-4000 

5 

5 

Horizontal 

1650 

2,750 

0 

-5000 

-5000 

-5000 

6 

10 

Horizontal 

4260 

7,100 

0 

-6500 

-6500 

-6500 

7 

10 

Hori zontal 

-4030 

-6,716 

0 

6500 

6500 

6500 

8 

10 

Horizontal 

-1232 

-2,054 

0 

6500 

6500 

6500 

9 

10 

Horizontal 

-3840 

-6,350 

0 

6000 

6000 

6000 

10 

10 

Horn zontal 

737 

1,229 

0 

6000 

6000 

6000 

11 

10 

Vertical 

3398 

5,663 

3000 

800 

-800 

-4000 

12 

5 

Vertical 

-143 

-238 

1800 

1400 

1200 

1200 

13 

20 

Horizontal 

-6692 

-11,154 

0 

6500 

6500 

6500 

14 

10 

Horizontal 

-3945 

-6,575 

0 

6500 

6500 

6500 

15 

10 

Horizontal 

-1740 

-2,900 

0 

-6500 

-6500 

-6500 

16 

10 

Horizontal 

2400 

4,000 

0 

-6500 

-6500 

-6500 

17 

5 

Vertical 

1901 

3,168 

3500 

2800 

2800 

4000 

18 

5 

Vertical 

-2770 

-4,616 

3500 

4800 

4800 

4000 

19 

5 

Vertical 

2640 

4,400 

2800 

2800 

2800 

4000 

20 

5 

Vertical 

814 

1,356 

1800 

1400 

1200 

1200 


a Load number 9.1 was a 1000°F (Mach 8) heating check. 

^Pressure load maintained at 0.7 5-1 b/i n2 . 

c Tests conducted at pressure loads of 0 and 0.75 1 b/ 1 n 2 . 

^Positive values indicate tension, negative values indicate compression. 
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Table 2 Test identification numbers and corresponding 
load conditions 


Test 

number 

Load 

condition 

Pressure, 

1 b/i n2 

Temperature, 

F° 

4.1.4 

4.1 

0.75 

70 

4.2.6 

4.2 

0.75 

70 

4.3.3 

4.3 

0.75 

70 

4.4.6 

4.4 

0.75 

70 

4.2.7 

4.2 

0 

70 

4.3.4 

4.3 

0 

70 

4.4.7 

4.4 

0 

70 

7.2.4 

4.2 

0.75 

550 

7.3.5 

4.3 

0.75 

550 

7.4.4 

4.4 

0.75 

550 

7.1.8 

7.13 

0.75 

550 

6.8.3 

3.8 

0.75 

550 

8.2.2 

4.2 

0 

1000 

8.2.2 

4.2 

0.75 

1000 

8.3.2 

4.3 

0 

1000 

8.3.2 

4.3 

0.75 

1000 

8.3.5 

4.3 

0 

1000 

8.3.5 

4.3 

0.75 

1000 

8.4.6 

4.4 

0 

1000 

8.4.6 

4.4 

0.75 

1000 

8.1.3 

8.13 

0 

1000 

8.1.3 

8. l a 

0.75 

1000 

Single-Panel 
Test 1.1 

N x * 0, N x y = 0 

0 

70 

Single-Panel 

N x * 0, N x y = 0 

0 

70 

Test 3.1 


. _ 

. 


a Equivalent to load condition 4.1. 
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Fig. 1 Geometry of tubular panel (dimensions in inches). Panel weight 
is 13 lb. 



E 29938 


Fig. 2 Photograph of Bene 41 tubular panel with 
heat shield supports removed. 
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170 x 10 3 



Fig. 3 Plasticity correction curves for local 
buckling of Rene 41 circular arc element, where 
T = 70°F, E = 31.6 x 10 6 lb/in 2 , and f oy = 

12 5 x 10 2 lb/m 2 (Ref. 2, eqs. 12-54 and 
12-55 setting n = 25). 


160 x 10 3 



Fig. 4 Plasticity correction curves for local 
buckling of Rene 41 circular arc element, where 
T = 550°F, E = 29.1 x ioS lb/m 2 , and f cy = 

120 x 10 s lb/m 2 (Ref. 2, eqs. 12-54 and 12-55 
setting n = 22.2). 


160 x 10 3 



Fig. 5 Plasticity correction curves for local 
buckling of Rene 41 circular arc element, where 
T = 1000° F, E = 26.5 x io 6 lb/m 2 , f ay = 

117 x io 3 lb/m 2 (Ref. 2, eqs. 12-54 and 12-55 
setting n = 18.5). 
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32 x 10 6 



Temperature, °F 


Fig. 6 Modulus of elasticity and modulus of rigidity 
as functions of temperature for Rene 41 (from Ref. 12). 



Fig. 7 Force/ stiff ness plot for local buckling. 
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E 26145 


Fig. 8 Hypersonic wing test structure (HWTS). 



DFRC80-236 


Fig. 9 Position of hypersonic wing test structure in proposed HR A. 
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Fig. 10 Hypersonic wing test structure with heat shields removed to 
show substructure . 


Up 



Fig. 11 Applied mechanical load distribution on HUTS and locations of five test 
tubular panels. 
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Fig. 12 HWTS combined mechanical and thermal loading teat setup, 



17 









E 40425-2 


Fig. 13 Tubular panel test assembly installed in the 
testing machine for axial compression buckling test. 
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Panel 1 


Panel 2 


Panel 3 


Wing root 


Panel 4 


Panel 5 


Fig. 14 
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a Rosette strain gages (RSQ) 

• Foil axial strain gages (ASG) 

■ Capacltatlve strain gages (CSQ) 

( ) Strain gages on inner (upper) surfaces 
I I Rosette strain gages used for elevated- 
temperature tests (panel center) 


Forward 

4 


Outboard 


Locations of strain gages on five tubular test panels (view looking downward). 
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E 33783 


Fig. 15 Four Rene 41 tubular panels attached to the hyper- 
sonic wing test structure for buckling tests ( view looking 
up and aft at the lower side of test structure) . Panel 3 
was removed to show the pressure pan. 
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